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Abstract: Many beyond the Standard Model (BSM) scenarios involve Higgs couplings to 
additional electroweak fields. It is well established that these new fields may modify /i — t- 77 
and h — )■ 7Z decays at one-loop. However, one unexplored aspect of such scenarios is that by 
electroweak symmetry one should also expect modifications to the hZZ coupling at one-loop 
and, more generally, modifications to Higgs production and decay channels beyond tree-level. 
In this paper we investigate the full BSM modified electroweak corrections to associated Higgs 
production at both the LHC and a future lepton collider in two simple SM extensions. From 
both inclusive and differential NLO associated production cross sections we find BSM-NLO 
corrections can be as large as ©(> 10%) when compared to the SM expectation, consistent 
with other precision electroweak measurements, even in scenarios where modifications to the 
Higgs diphoton rate are not significant. At the LHC such corrections are comparable to 
the involved QCD uncertainties. At a lepton collider the Higgs associated production cross 
section can be measured to high accuracy (0(1%) independent of uncertainties in total width 
and other couplings), and such a deviation could be easily observed even if the new states 
remain beyond kinematic reach. This should be compared to the expected accuracy for a 
model-independent determination of the /i77 coupling at a lepton collider, which is 0(15%). 
This work demonstrates that precision measurements of the Higgs associated production cross 
section constitute a powerful probe of modified Higgs sectors and will be valuable for indirectly 
exploring BSM scenarios. 
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1 Introduction 

Last year the ATLAS and CMS collaborations discovered a new scalar in the bosonic Higgs 
search channels [1, 2] exhibiting properties consistent with the Standard Model (SM) Higgs 
boson with mass nih — 125 GeV. This mass is fortuitous from both experimental and theo- 
retical perspectives since multiple Higgs decay channels can potentially be observed, allowing 
for a multi-facetted exploration of electroweak symmetry breaking. The rich experimental 
possibilities, combined with the fact that the Higgs is the focus of many new physics scenar- 
ios, propels the Higgs to center stage in the effort to connect BSM theory with experiment. 
This leads to a number of phenomenologically important questions. If there is new physics, 
what is its nature?^ How might we hope to observe it directly or indirectly? Addressing the 
latter question requires understanding in detail how new physics scenarios might affect Higgs 
production and decay properties and, consequently, precision analyses at present and future 
colliders. Even if no significant deviations in the Higgs properties are observed from the first 
run of the LHC, these issues remain at the top of the phenomenological agenda. 

^Due to early tentative hints for an excess in /i — >■ 77 rates this question has stimulated a great deal 
of theoretical activity. See e.g. [3, 4] for explicit models and also discussions of the implications of BSM 
contributions to the /177 coupling. These hints have persisted in recent ATLAS analyses [5], but only at the 
level of 2a, and CMS have not yet released an updated analysis [6], so the status of this excess is very uncertain 
and we do not consider it a motivating factor in this work. 
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Figure 1: New electroweak-charged fields coupled to the Higgs contribute at one loop to 
(a) /i — )• 77 decays and (b) h — )• 7Z decays. The branching ratios to these final states are 
sensitive to the total Higgs width, which depends on Higgs couplings to other SM and BSM 
fields. The total rate in these channels also depends on the Higgs production cross section. 
Thus experimental determination of the h'yy and h'jZ couplings, either at the LHC or a 
lepton collider, is subject to uncertainties in all of the Higgs couplings. For these reasons, 
at a 250 GeV lepton collider with 250 fb~^ integrated luminosity the /i77 vertex can be 
determined to an accuracy of 0(15%) [7]. 
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Figure 2: Higgs associated production. New electroweak-charged fields coupled to the Higgs 
typically contribute to this amplitude at one loop, interfering with the tree-level amplitude. 
The hZZ coupling can be measured to a high degree of accuracy (0(1%) [7]) at a lepton 
collider by measuring Z-recoils in the inclusive associated production process, regardless of 
uncertainties regarding the total Higgs width. Hence a lepton collider gains better sensitivity 
to the hZZ coupling than to the /i77 and hjZ couplings, and opens up the possibility to 
discover new physics through anomalous contributions to the hZZ vertex. 



Perturbative BSM scenarios which involve Higgs couplings to new electroweak fields are 
often characterized by new states charged under SU(2)i x U(l)y which obtain, at least par- 
tially, their mass from the Higgs sector. In Figs, la, lb, and 2, we schematically demonstrate 
that by electroweak symmetry Higgs /i — )• 77 decays, h — )• 7Z decays, and associated pro- 
duction at the Tevatron, LHC, and a future lepton collider^ are typically influenced by such 

^We will use the term 'lepton collider' in reference to any future high energy e'^e^ synchrotron, linear 
collider or Higgs factory. 
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couplings. In particular, new physics contributions that enter h ^ and h — )• Z7 at leading 
order are directly correlated with the next-to-leading order (NLO) electroweak corrections to 
Higgs boson production in associated production.'^ This observation is not merely a technical 
curiosity, but is relevant across the broad scope of the experimental Higgs frontier since the 
associated production process lies at the heart of a number of critical measurements, which 
we now briefly review.^ 

While there are emerging hints of /i —>• rr at the LHC [10, 11] and /i — )■ 66 at the 
Tevatron [12], the Higgs boson candidate remains to be discovered in the fermionic final 
states. Given that the decay /i — )■ 66 in the SM for rrih — 125 GeV is the largest decay 
channel, with a branching ratio BKsAiih — t- 66) ~ 60%, not observing the Higgs boson in 
this channel indirectly induces a large uncertainty on the extraction of Higgs couplings to 
other fields. We note however that this uncertainty can be artificially decreased by making 
assumptions on the total width of the particle. Looking to the future, observing /i — t- 66 at the 
LHC in either gluon fusion or weak boson fusion will be difficult, due to overwhelmingly large 
QCD backgrounds, trigger issues, or the involved selection criteria which in the latter search 
channel typically involve central jet vetoes [1, 2, 13, 14].^ However, given the success of jet 
substructure techniques on boosted final states from associated Higgs production pp — )• hZ 
[15, 16], there is good reason to believe that such an analysis will shed light on — )• 66 in the 
future [17-19]. Thus good theoretical control over the associated production cross section is 
integral to future determination of the fermionic Higgs couplings at the LHC.^ 

The status of the associated production process is elevated further at any future lepton 
collider since at low energies (~ 250 GeV) associated production is the dominant Higgs 
production process and influences all Higgs observables. Furthermore, measurements of the 
inclusive associated production cross section will allow determination of the hZZ coupling 
to 0(1%) accuracy, regardless of the total Higgs width, making this the most accurately 
determined of all the Higgs couplings. The associated production cross section also has 
ramifications for precision determination of a number of Higgs couplings as the total Higgs 
decay width can be inferred through a combination of the inclusive associated production 
cross section, specifically Z — )• ix^ ^~ recoiling against the Higgs, and other production and 
decay channels. Hence, if there are any consequences of a modified Higgs sector for associated 
production, they will be important at lepton colliders. 

Non-decoupling parts of the one-loop contributions of Figs, la and lb are sensitive to 
physics beyond the SM even if the new mass scale lies well above the LHC or lepton collider 

^New electroweak fields will also contribute in the Z self-energy and 7/Z mixing. We include these contri- 
butions in all calculations, but not in diagrams. Similar corrections will also enter into h — >■ ZZ* decays. We 
do not consider this here, although it is certainly deserving of further investigation. 

*For other BSM applications of the associated Higgs production process see e.g. [8, 9]. 

^Although difficult, we note that observation of /i — > 66 in the weak boson fusion channel is possible in the 
future. We thank Markus Klute for bringing this to our attention. 

^It should also be noted that one may be able to study different production cross sections with the same 
Higgs decays. For example, h — ^ WW* decays may be observable separately in gluon fusion, vector boson 
fusion, and associated production channels (see e.g. [20, 21]) leading to relations between cross sections. 
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energy reach. If the new physics mass scale A is large enough, it is customary to parametrize 
deviations from the SM by including higher dimensional operators that arise from integrating 
out the dynamics above A in an effective field theory language.'' This approach has the 
benefit of being model independent, but is limited due to the very nature of effective field 
theory. This does not only follow from the non-renormalizable character of effective field 
theory, which systematically hinders theoretical precision [24], but sheer experimental reality 
puts the link between experiment and theoretical interpretation under stress. Many well- 
motivated physics extensions, especially in the electroweak sector, result in a phenomenology 
that can be consistent with current data without having a too high new physics scale A ~ 
100 GeV. In this case, a serious drawback of the effective theory approach is that experimental 
signal versus background discriminating selection criteria typically need to be designed to 
probe momentum transfers ~ A simply in order to overcome the huge SM backgrounds. 
This means we would probe the effective theory at scales at which it cannot be considered 
valid. An example for such an analysis is again the boosted /i — t- 66 associated production 
where pt{Z) > 120 GeV facilitates both triggering the event and rejecting the large tt+jets 
background [15, 16]. Therefore the precision analysis of SM extensions at similar energy 
scales can become model-dependent and only dedicated calculations allow one to draw a 
qualitatively complete picture.® 

For the aforementioned reasons we choose to study particular models, since although 
one loses some degree of generality, this enables precision calculations and also serves to 
illustrate the features one might expect in a broader range of SM extensions. To determine 
where the new electroweak physics may become apparent, in this paper we investigate the 
correlation between two obvious observables in which new electroweak physics may lead to 
modifications at one-loop: corrections to /i — ?■ 77 decays, which is very commonly studied, and 
NLO electroweak corrections to associated Higgs production at the LHC and a future lepton 
collider where this cross section can be measured to high precision. To do this we employ 
two concrete models. We first discuss a vector- like leptonic fourth generation in Sec. 2.1 and 
then a charged scalar Higgs sector extension of the SM in Sec. 2.2. For both models we 
establish consistency with current constraints in terms of the precision electroweak S, T, and 
U, parameters [25, 26]. 

Simple electroweak extensions typically involve singlets under Quantum Chromodynam- 
ics (QCD) SU(3)c- Depending on its spin, color-charged matter can compete with SM elec- 
troweak fields at one-loop, but can also be directly constrained at the LHC due to limits 
on the production of colored particles [27]. Discarding non-trivial color assignments, which 
we will do in the following, does result in some degree of qualitative theoretical deficiency, 
however if one wished to include colored matter then, by considering Fig. la, Fig. lb, and 
Fig. 2, it is clear that the BSM amplitude simply scales uniformly with the casimir of the 
SU(3)c matter representation. 

'^For early work in the gauge and Higgs sectors see Refs. [22, 23] 

**Note also that by typically probing a large center-of-mass energy the information encoded in electroweak 
precision observables such as the 5", T, U observables [25, 26] do not govern the collider phenomenology. 
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In Sec. 3 we consider electroweak corrections at lepton colliders, at which the impact 
of the electroweak quantum effects is most pronounced and their observation is challenged 
least. We set the stage in the context of the SM in Sec. 3.1, bridging our work to the existing 
literature. In Sec. 3.2 we discuss modifications of the NLO electroweak phenomenology in the 
presence of BSM interactions. For completeness we provide a survey of the expected BSM 
NLO electroweak effects at the LHC in Sec. 4, utilizing exemplary parameter points that are 
motivated from our lepton-collider discussion of Sec. 3.2. We summarize this work and give 
our conclusions in Sec. 5. 

2 Example BSM Scenarios 

2.1 Vector-Like Fourth Generation Leptons 

As a first illustrative BSM scenario we discuss a model with an additional family of vector-like 
leptons.^ The model introduces a full vector-like generation of leptons and the Lagrangian, 
suppressing the kinetic and gauge interactions, is^° 

— £ D mi^i'j^i'^ + mee[e^ + m^u'li^'p^ + h.c. (2.1a) 
+ YMH)e'R + Y:,{FLia'H^y^ + Y^{e%H)e'i + Y^ii'^ia^H^yi + h.c. (2.1b) 

where f^,f^ = (2,-1/2), e'[,e'j^ = (1,-1), and z^^,z^^ = (1,0) under SU(2)l x U(l)y. The 
left-right symmetric character of this extension guarantees cancellation of anomalies. 

The parameter space of this model is large, and so we choose to make some restrictions to 
facilitate exploration of the relevant phenomenology. In the following we assume all vector-like 
masses take a common value 

rrii = rUe = rriy = my , (2.2a) 

where the subscript V simply denotes that this mass contribution is vector-like. We also 
assume that all charged fields have the same Yukawa couplings and obtain the same mass 
contributions from the Higgs. We make a similar simplification for the neutral fields but allow 
their Yukawa couplings to differ from the charged fields. Thus we define 

Y^vh/^j2 = Y^vh/V2 = mch (2.2b) 
Y^Vh/V2 = Y^Vh/V2 = mch + ^u , (2.2c) 

where the subscript Ch denotes that this mass contribution is chiral in nature. This simplified 
version of a vector-like fourth generation of leptons is hence fully specified by the parameter 
set {mch,mv,A^}. 

^In Ref. [3] the model has been discussed in the context of enhanced /i — >■ 77 decays, dark matter and 
constraints from direct searches. 

^"PMNS-like mixing of the leptons with the first three generations is quite constrained and can be avoided 
by imposing global symmetries. This is implicitly assumed in the following. We borrow the conventions of 
Ref. [3]. 
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Finally we define to be the mass of the lightest charged fermion and throughout 
will assume niEi > GeV, just beyond the direct production reach of a 250 GeV lepton 
collider. 

The mass eigenstate basis that results from diagonalizing (2.1) with (bi)unitary transfor- 
mations generates mixings in the gauge and Higgs interactions (2.1b). While these mixings 
are important for processes involving production and decay of the new fields, they will also 
typically leave footprints in the low energy physics, below the energy scales at which the 
new physics is directly accessible. While the focus in this work is on modifications of Higgs 
physics at scales A ~ mh, we must also consider the effects on physics at scales A ~ mz- Such 
modifications are conventionally expressed in terms of the precision electroweak S, T, and U, 
parameters [25, 26], which are sensitive to the degree of the model's non-chirality and any ad- 
ditional custodial isospin violation. For our parameter identifications Eqs. (2.2a)-(2.2c), the 
former is controlled by mcfi and the latter by Aj^. In the language of effective field theory, 
U follows from a dimension eight operator as opposed to the dimension six operators which 
parametrize S and T, and hence is typically suppressed [28]. In this model contributions 
to S and T can be considerable, and so when considering the influence of the new fields on 
Higgs associated production and decays we will also include precision electroweak constraints 
throughout. 

2.2 New Electroweak Scalars 

As a second illustrative BSM scenario we consider a model containing additional scalars. In 
particular we include an extra scalar doublet, (j), with the same quantum numbers as the Higgs 
doublet. For the sake of simplicity we assume (p is charged under an approximate global U(l) 
symmetry^^ which restricts the scalar potential to contain the following terms 

V D ml\cl)\'' + X\H\^\<P\^ + X'\H ■ <p^^ . (2.3) 

After electroweak symmetry breaking the Higgs vev leads to a mass splitting between the 
neutral and charged scalars, (po and (j)-^ and also generates trilinear interactions between the 
Higgs and pairs of scalars. It is useful to trade the parameters of the above Lagrangian, 
{m,^, A, A'}, for the more intuitive set {m<^^, A(j,_^, A^}, where is the mass of the charged 
scalar, A^^ is the trilinear coupling between the Higgs and charged scalar C D Atp^h\4>^\'^ , 
and is the mass splitting between the neutral and charged scalars 

A<^ = - m^^ . (2.4) 

^'^Due to the sterile neutrino extension, the model also allows additional Majorana mass terms, which can 
lead to a good dark matter candidate. For the present context, i.e. the ILC and LHC collider phenomenology 
at NLO, the only impact of the Majorana masses is a modified mixing of the gauge interactions, which has no 
effect on the observed /i — )■ 77, and we do not consider Majorana mass terms in the following. 

^^We assume this global symmetry is broken by some small amount to allow for </> to decay, with the proviso 
that the operator responsible for this is small enough to be irrelevant for the loop corrections we discuss. 
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Figure 3: SM electroweak one-loop corrections to the associated production process at a 
lepton collider, (a) Total associated production cross section at a linear collider as a function 
of ^/s and (b) differential cross section at y/s = 250 GeV. In both panels the tree-level result 
is plotted in red (middle of a set of three), the result with only weak corrections included in 
blue (top of three), and the full one- loop result in black (bottom of three). The various Z- 
boson polarizations are plotted in continuous (unpolarized), dashed (longitudinal) and dotted 
(transverse) lines. 



As for the vector-like lepton model we will also include precision electroweak constraints on 
this model throughout and assume m^^ > 125 GeV, such that the lightest charged scalar is 
beyond direct production reach at a 250 GeV lepton collider. 

3 (B)SM Electroweak Radiative Corrections at a Lepton Collider 

Although the Higgs was only discovered recently, the program to understand and calculate 
radiative corrections in various Higgs production and decay processes was initiated long ago. 
Measuring electroweak radiative corrections typically requires a precision machine, so in this 
section we focus on Higgs production at a lepton collider, and in particular consider an 
e~^e~ collider operating at y/s = 250 GeV, for which the dominant production mechanism is 
associated production, e'^e~ — )• hZ [8, 29]. Radiative corrections to associated production 
at a lepton collider were calculated some time ago [30-33], and before proceeding to BSM 
scenarios we revisit these calculations armed with known values for the top quark and Higgs 
mass. 
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3.1 SM 



Electroweak radiative corrections to e^e 



hZ (and also qq — )• hZ) consist of self energy 



corrections, vertex corrections and box diagrams. We follow [33] and calculate these correc- 
tions using the FeynArts, FormCalc, and LoopTools packages [34, 35]. We calculate 
counter-terms with dimensional regularization in the complete on-mass-shell renormalization 
scheme [36-39] and also split the radiative corrections into two classes, electromagnetic and 
weak, denoting the full associated production cross section 



where ctq is tree-level result, aweak contains the weak corrections which do not involve photons 
within loops, and (Jem includes the virtual photon correction to the eeZ vertex as well as 
soft photon emission. Including a massless photon in the virtual contributions of ctem we 
need to be inclusive on the final state to guarantee the cancellation of infrared singularities 
according to the Kinoshita-Lee-Nauenberg theorem [40, 41]. In the soft photon limit the real 
photon emission (bremsstrahlung) contribution factorizes to the born cross section [42] and we 
include the corresponding factor in our definition of uem for photon energies E{^) < 0.1^/s, 
again following [33] . Hard photon corrections at the given order of the perturbative expansion 
are insensitive to new physics contributions and we do not consider them in the following; 
they can be vetoed in the actual analysis. To get a feeling for the magnitude of contributions 
coming from loops of heavy particles we can also investigate the weak corrections alone by 
defining 



with the electromagnetic corrections omitted. We have verified the cancellation of UV di- 
vergences for both the weak and electromagnetic corrections and also the cancellation of IR 
divergences in the electromagnetic corrections. 

Comparing our results with those presented in [33] we find excellent agreement. Updating 
with the known top-quark and Higgs mass (which we set to nih = 125 GeV throughout) we 
calculate the weak corrections within the SM at a future 250 GeV lepton collider. In Fig. 3a 
we plot the total associated production cross section as a function of y/s (which peaks near 
^/s « 250 GeV) and in Fig. 3b we plot the differential cross section at ^/~s = 250 GeV. The 
weak corrections are typically small and positive for all polarizations and scattering angles. 
The electromagnetic corrections are slightly larger and negative, thus, when combined, the 
weak and electromagnetic corrections lead to an overall reduction in the associated production 
cross section. The weak corrections remain below 5% of the tree-level cross section, and for 
unpolarized Z-bosons constitute an 0(3%) correction at ^/s = 250 GeV. Although these 
corrections are small it is expected that the hZZ coupling can be determined to 0{1%) 
accuracy at a 250 GeV lepton collider with 250 fb~^ integrated luminosity [7], and so weak 
radiative corrections can in principle be measured at such a machine. 



f = (To + aEM + Crweak 




aw = 0-0 + (TWeak 
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It is clear that SM weak radiative corrections, although small, are large enough to be 
relevant to lepton collider Higgs studies. Their magnitude, 0(3%), would lead one to expect 
corrections arising due to perturbative BSM scenarios to be similar. 

3.2 Beyond the SM 

We calculate these corrections following the same methods as for the SM weak corrections 
and include all contributing diagrams, including those that contribute to vector boson self 
energies. We define the BSM correction to the associated production cross section 



crBSM{e~^e Zh) - asM{e~^e Zh) 



o■o(e^ 



Zh) 



and in a similar manner define 



'h-yy 



'^BSMjh 77) _ BRb5m(/^ 77) 
^SM{h 77) BRsM{h 77) 



ptot 
BSM 
ptot 
^ SM 



(3.3) 



(3.4) 



Recent ATLAS fits indicate total width constraints T^^j^^ 



/r5M ^ 2.12 [43]. While this 
number results from an overall scaling factor of LHC measurements and will therefore be 
limited systematically, the total Higgs width can be accessed at a 250 GeV lepton collider 
via a measurement of the e^e~ — )• I'uh channel with subsequent decay h — )• W~^W~ alone at 
the 10% level [44]. We will assume throughout that all new fields are too heavy to introduce 
important new Higgs decay channels, and hence ^^bsm ~ ^^SM^ such that modifications of 
the diphoton decay rate arise solely due to BSM contributions to the decay amplitude. 

In general BSM contributions to 6azh and Rh^yy are correlated, however the form of this 
correlation is model dependent; the dependence arising through the spins and electroweak 
charges of the additional BSM fields. For this reason we consider this correlation in the two 
sample BSM scenarios described in Sec. 2. 

3.2.1 Vector-Like Fourth Generation Leptons 

This simplified model, presented in Sec. 2.1, is defined by the parameters {mch, 'mv, ^u}- In 
each panel of Fig. 4 we fix A^, to a specific value and vary mch and my , plotting contours of 
the lightest charged fermion mass, niE^, the Higgs diphoton decay rate Rh-yy, and corrections 
to the associated production cross section 5a zh- We also shade regions disfavored by precision 
electroweak measurements. 

Since modifications of the Higgs diphoton decays do not involve the neutral fields, Rh"/-y 
depends only on mch and my- It is clear that over the majority of allowed parameter space 
corrections to Rh-yy are typically quite modest, C'(10's%) compared to the 0(15%) model- 
independent experimental accuracy of the /177 coupling. Achieving an enhancement of the 
diphoton decay rate requires mv > mch- This might appear surprising since this implies the 
dominant mass contribution does not come from the Higgs. However, if the fermion mass 
comes solely from the Higgs the diphoton decay rate never receives an enhancement since the 
contributions from loops of fermions are always smaller in magnitude and opposite in sign 
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Figure 4: Contours of the modified Higgs diphoton decay rate, Rh-yy (dashed red), and 
BSM corrections to the associated production cross section dazh (sohd black) as defined in 
Sec. 3.2, for the vector-like fourth lepton generation model. Contours of the Hghtest charged 
fermion mass, mcj, are also shown (dotdashed purple) but not labelled. For our restricted 
parameter choice we have niEi = \my — mch] so the lightest fermion mass can be read off 
from the intersection of these lines with the x or y axes, (a) corresponds to the case where 
the Higgs contributions to charged and neutral fermion masses are equal, Ai, = GeV, (see 
Eq. (2.2c)) and (b) to the case where the neutral fields couple more strongly to the Higgs 
than charged fermions, Ap = 50 GeV. Regions where precision electroweak corrections exceed 
measured values by la are shaded, with A^* in gray and AT in blue. Parameter values where 
modifications to the diphoton rate are not significant {e.g. 0{'^ 2a) corresponding to the 
region to the left of the Rh'yy = 1.3 contour) typically imply a reduction in the associated 
production which can comfortably exceed 1% compared to the SM prediction. This would 
be observable at a 250 GeV lepton collider which can measure this cross section to 0(1%) 
accuracy. 

to the dominant amplitude from loops of VF-bosons.^'^ If one increases my then the lightest 
fermion can be kept light by also increasing rrich- This then allows for much larger couplings 
to the Higgs for a given light fermion mass and hence a larger contribution to the Higgs 
diphoton amplitude. 

Considering Fig. 4 it is clear that within this model BSM-NLO corrections typically lead 
to a reduction in the associated production cross section of 1% or greater, even in cases 

^^To achieve a large enhancement in this scenario one could always include multiple copies or larger charges, 
however such modifications are not within the definition of the model being studied here. 
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[GeV] 

Figure 5: Precision electroweak constraints on the electroweak scalar model. By far the 
strongest constraint comes from the T-parameter and consistency at la with precision elec- 
troweak measurements requires [A^l < 84 GeV. Inconsistency with the 5-parameter only 
arises in small wedges of parameter space. We also require that the neutral scalar has mass 
greater than m^g > 110 GeV to ensure consistency with LEP bounds. 



where corrections to the diphoton rate are not significant. With sensitivity to the associated 
production cross section at 0(1%) such deviations would be easily measurable at a lepton 
collider, even if the Higgs diphoton rate is within 0(1 — 2)a of the SM value. 

Although not shown here we have also investigated the dependence of the cross section 
corrections on Z-boson polarizations and scattering angles. The corrections to the associated 
production cross section do not discriminate much between longitudinal and transverse Z- 
boson polarizations, with similar corrections to both. Furthermore, the corrections are also 
largely independent of the scattering angle. 

In conclusion, if the Higgs is coupled to additional heavy charged fermions then, even 
if these fermions remain beyond direct collider reach, and even if their corrections to the 
diphoton decay rate are not significant, their presence should be evident at a lepton collider 
through modifications to the associated production cross section. 

3.2.2 New Electroweak Scalars 

This model is defined by the parameters {m^_^_ , A^j,^ , A^} and for a particular choice of these 
parameters both 5azh and -R/177 are completely determined. As shown in Fig. 5 the charged- 
neutral scalar mass splitting A,^ is constrained by precision electroweak measurements to 
be lA^I < 84 GeV. With this constraint imposed this model is consistent with precision 
electroweak measurements for the full range of rrifp^. 

In Fig. 6 we plot contours of 6azh and Rh^-y for fixed values of A,^ and varying values 
of {m^^,A^^}. From these plots it is clear that, within this particular model, it is possible 
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Figure 6: Contours of 6azh (see Eq. (3.3)) in solid black, and Rh^yy (see Eq. (3.4)) in dotted 
red, for the charged scalar model with varying scalar masses and Higgs couplings. Corrections 
to the associated production cross section are often 0(few)% even when corrections to the 
diphoton rate are not significant {Rh^yy — 1| ^ 0.3, (corresponding to an 0{2a) deviation in 
Rhy-y)- With sensitivity to the associated production cross section at the 0(1%) level these 
modifications to associate production could be observed at a linear collider. 

to generate deviations in the associated production cross section which are larger than the 
0{1%) achievable experimental resolution. Again, although not shown here, the corrections 
to the associated production cross section do not discriminate greatly between longitudinal 
and transverse Z-boson polarizations and are also largely independent of the scattering angle. 

Hence, similarly to the model with vector-like fermions, if the Higgs is coupled to new 
charged scalars close to the weak scale then, even in the case where these scalars remain 
beyond direct reach of the LHC and lepton colliders, there is still the potential for measuring 
significant deviations in the associated production cross section at a linear collider, assuming 
the cross section can be measured to an accuracy of 1%. 

3.2.3 Dominant Corrections in the Scalar Model 

In Fig. 6 one can see that the one- loop BSM corrections, 6a zh, are, to a reasonable approxi- 
mation, a quadratic function of the trilinear coupling Afp_^ . Of all contributing diagrams, only 
the scalar loops which contribute to the Higgs wavefunction renormalization have a quadratic 
dependence on the scalar coupling. This wavefunction renormalization then feeds into the 
hZZ vertex through a field redefinition to give canonically normalized external states. This 
contribution to the effective vertex is finite and if we only include it then modifications to the 
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associated production cross section can be simply written as 

where t_|_ = {mh/2m^^)'^ ^-^"^ ''"o = ('^/i/2m0Q)^, which are commonly found quantities in loop 
integrals. A^^ is the trilinear coupling between the Higgs and charged scalars, as defined in 
Sec. 2.2, and is the trilinear coupling between the Higgs and neutral scalars, which can 
also be written as ^ 

^00 = ^4>+ + —("^00 - '^1+) ' (3-6) 
where Vh = 246 GeV is the usual Higgs vev. Finally we define the function 

^(^) = 1 + ^7^^=iy 0°^ (l - 2t - 2 Vr(r - 1)) - log (l - 2r + 2^t{t - 1)) ) , (3.7) 

which is related to the derivative of the scalar two-point function. Since < r+ o < 1 for the 
cases of interest the denominator and arguments of the logarithms in Eq. (3.7) are complex, 
however the final result is real. 

Comparing with the full one-loop corrections we find that whenever the scalar coupling 
is large, A^^^ > Vh, then Eq. (3.5) gives a reasonably good approximation to the full ex- 
pression. It is interesting that the dominant one-loop correction is entirely independent of 
the electroweak charges of the scalars, demonstrating that the associated production cross 
section is a sensitive collider-probe of scenarios where the Higgs is coupled to neutral scalars 
via the Higgs-portal [45-52]. This is especially interesting as it applies even if the scalars 
are too heavy to introduce invisible decays. Thus, if there are no modifications to the Higgs 
diphoton decay rate, and if the Higgs has no invisible decay width, the associated production 
cross section can be used to probe the influence of additional neutral states at one loop.^^ 

4 BSM Electroweak Radiative Corrections at the LHC 

We now turn to a discussion of the extended electroweak corrections at the LHC. Obviously, 
the electroweak NLO corrections [53] are not as important as the involved higher order and re- 
summed QCD corrections [53-59] which currently give rise to an 0(10%) residual perturbative 
uncertainty [60]. It is however not unreasonable to expect that this number will improve over 
time when more data becomes available and theoretical development is stimulated, such that 
potentially large electroweak modifications can be resolved. 

Technically, we proceed along the lines of Sec. 3 to calculate the electroweak NLO cor- 
rections. We use the CTEQ611 parton distribution functions (pdfs) [61] following Ref. [53] 
within the Vbfnlo framework [62] to handle the pdf-folding and phase space integration for 
^-£-(7) < 0.1 ^f§ {s denotes the partonic s). We have again checked the cancellations of 

^''As this leads to an almost universal correction to all Higgs couplings evidence for such a scenario might 
also arise through studying other Higgs production channels. 
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Figure 7: (a) transverse momentum, and (b) rapidity, differential distributions of the Z 
boson in pp — )• hZ production at the LHC y/s = 14 TeV. We show leading order (LO) 
distributions as well as electroweak corrections in the SM and the BSM models discussed in 
the previous sections. The BSM parameter choices are given in Eq. (4.1) and Eq. (4.3). The 
integrated cross section receives corrections Sazh ~ —15% in the charged scalar model and 
^^Zh — —8% in the vector-like fermion model. The correction relative to the SM value, as a 
function of pT, is shown separately in Fig. 8. 



both UV and IR singularities and we find excellent agreement comparing to integrated re- 
sults from FormCalc/LoopTools which validates our phase space integration. Our leading 
order matrix element has been cross checked against MadGraph [63] for individual phase 
space points. 

As done for the e^e~ case we regulate the soft photon divergencies with a slicing pa- 
rameter which is equivalent to dimensional regularization [64], and we do not consider either 
initial state photon contributions or hard photon corrections.^^ 

In Fig. 7 we show differential distributions for the SM and our BSM scenarios. For a 2 — t- 2 
scattering process the most relevant distributions are transverse momentum pT and rapidity 
y. Dependence on the azimuthal angle, is trivial. The remaining collider observables follow 
as functions of {pT,y,^) and the involved masses mz,Tnh. 

^^Electroweak processes also typically depend on the electroweak scheme that is used to construct the 
parameters of the electroweak sector from a minimal set of observables [53]. This dependence can be large but 
is not important for us since it does not involve the BSM contributions. 

^^In principle, the presence of coUinear mass singularities requires extending the list of subprocesses to 
photon-induced processes and the DGLAP kernels accordingly. Such pdf sets exist, see e.g. Ref. [65], but their 
contribution to the present process is known to be negligible < 1% for our purposes [66, 67]. Note that none 
of these contributions is sensitive to the BSM extension that we study. The full treatment will not change our 
results; especially the adoption of Eq. (3.3) to the LHC situation is insensitive to all non-BSM effects. 
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The charged scalar mass scenarios we will study in the following are given by 



(a) {m^^ A^^, A^) = (250 GeV, -2.0 TeV, 80 GeV) , (4.1) 
(a) (m<^+ A^^, A^) = (250 GeV, +0.7 TeV, 80 GeV) , (4.2) 

and we choose a parameter point for the vector-like lepton scenario 

(b) (mcfc, my, A^) = (500 GeV, 1000 GeV, GeV) . (4.3) 

These parameter points are motivated from our results of Sec. 3.2, Figs. 4 and 6. We choose 
(a) and (b) in the lepton collider {Rh-y-y, a{hZ)} correlation such that the diphoton decay rate 
takes an extreme value Rhj^y — 1-5, while (a') is a parameter point with a ~ 20% decrease in 
Rh'yy compared to the SM.^^ 

We find LHC BSM cross sections 
6a{hZ) ~ —15%, —8% for the parameter and 
model choices (a) and (b), respectively. For 
(a') we obtain 6a{hZ) ~ —1%. These numbers 
qualitatively reproduce the lepton collider re- 
sults of Sec. 3.2 (we remind the reader that 6a 
is mostly driven by the BSM corrections). As 
in the lepton collider case, the scalar extension 
shows a larger deviation as compared to the 
fermion model. The scalar model is also less 
constrained by electroweak precision results, 
which follows from a formal decoupling of the 
non-oblique corrections (see also Sec. 3.2.3). 

In Fig. 8 we show the phase space depen- 
dence in terms of the Z boson transverse mo- 
mentum, which is sensitive to the BSM con- 
tributions, especially when we start probing 
the new physics mass scale. This justifies our 
choice to work in concrete models as compared to effective theories mentioned in the intro- 
duction. This effect is more clearly visible in Fig. 8. The total cross section is dominated 
from low pt configurations, yet for analyses where associated production is accompanied with 
a high pt cut, such as in /i — )• bb where pt{Z) — 150 GeV, we observe a turn-on of the new 
physics contribution which further decreases the cross section compared to the SM expecta- 
tion on the inclusive level. The phase space dependence of the correction, however, is highly 
dependent on the model. The overall effect is small and in the percent range, but should 
nonetheless be considered in precision studies that aim to determine Higgs properties at the 

^'^The fourth generation vector-like lepton model does not admit Ru-y-i < 1 due to electroweak precision 
constraints 
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BSM electroweak corrections of the con- 
sidered models and parameter points as 
a function of the Z boson transverse mo- 
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LHC. Within the context of the charged scalar model the parameter point (a') does not result 
in important effects on LHC associated production. 

We have only discussed pp — )• hZ but expect our results should generalize to pp — )■ hW^ 
qualitatively as in the SM [53], and so these results also have relevance for the pp — )• hW — )• 
WWW channel [20]. 

5 Conclusions 

In this work we have investigated the higher order electroweak corrections to associated Higgs 
production at a future lepton collider and the LHC in electroweak extensions of the Higgs 
sector. Because higher order corrections are model and phase space dependent, we have 
studied them within two simple BSM scenarios which involve vector-like leptons and charged 
scalars. Our results are of importance for the consideration of Higgs physics at a future lepton 
collider since couplings and branching ratios can be measured to high precision. Reconciling 
the new physics contributions with current limits from LEP measurements we find regions of 
parameter space with 0(10%) or larger deviations of the Higgs associated production cross 
section, in comparison to the SM expectation. These corrections are much larger than the 
expected lepton collider precision, and are possible even when corrections to the diphoton 
rate are within 0(1 — 2)a of the expected precision. 

We find similar corrections for LHC collisions. The phase space dependence is typically 
such that the integrated BSM cross section modifications are dominated by soft events, but 
high pt events exhibit the largest deviations. These corrections, when compared to the SM, 
can again be 0(10%), which might be resolvable for high luminosity studies. The phase space 
dependence of the BSM contribution then becomes a non-negligible parameter in precision 
cross section and coupling extractions, which, depending on the data, might challenge an 
effective field theory based analysis as de facto considered in Ref . [68] . 

While signals for any BSM physics at the LHC remain elusive it is increasingly important 
to explore the full breadth of potential experimental BSM signatures. Although the direct 
observation of new BSM states still remains a possibility, it may be that indirect evidence in 
the form of anomalous properties of SM states and processes might also uncover new physics. 
Since many well-motivated BSM scenarios are concerned with the Higgs sector, the newly 
discovered Higgs then becomes, perhaps, the strongest candidate for indirect hints of new 
physics and it is timely to understand the full complement of possible BSM modifications 
to its properties. This effort is already underway with intense study in various respects, 
however this work highlights that as we progress into the Higgs precision era it is important 
that theoretical progress in SM precision Higgs calculations comes accompanied by BSM 
precision Higgs calculations. The scope of BSM precision Higgs physics is extensive and in 
this work we have only considered corrections to Higgs associated production in two simplified 
models. However, this has demonstrated that if the Higgs is coupled to new electroweak states 
then there exist regions of parameter space where modifications to /i — t- 77 rates may remain 



- 16 - 



below statistical significance, while corrections to associated production, particularly at a 
lepton collider, may in fact lead to the most significant evidence for a modified Higgs sector. 
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